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Abstract 

A quantum absorption refrigerator driven by noise is studied with the purpose of determining 
the limitations of cooling to absolute zero. The model consists of a working medium coupled 
simultaneously to hot, cold and noise baths. Explicit expressions for the cooling power are obtained 
for Gaussian and Poisson white noise. The quantum model is consistent with the first and second 
laws of thermodynamics. The third law is quantified, the cooling power J c vanishes as J c oc T°, 
when T c — > 0, where a = d + 1 for dissipation by emission and absorption of quanta described by 
a linear coupling to a thermal bosonic field, where d is the dimension of the bath. 

PACS numbers: 03.65.Yz,05.70.Ln, 07.20.Pe,05.30.-d 



1 



I. INTRODUCTION 



The adsorption chiller is a refrigerator which employs a heat source to replace mechanical 
work for driving a heat pump [1]. The first device was developed in 1850 by the Carre 
brothers which became the first useful refrigerator. In 1926 Einstein and Szilard invented an 
abso.pt.on refrigerator with no .aoviag parts fl. This idea has heen incorporated recent* 
to an autonomous quantum absorption refrigerator with no external intervention {2, 4| . The 
present study is devoted to a quantum absorption refrigerators driven by noise. The objective 
is to study the scaling of the optimal cooling power when the absolute zero temperature is 
approached. 

This study is embedded in the field of Quantum thermodynamics, the study of thermody- 
namical processes within the context of quantum dynamics. Historically, consistence with 
thermodynamics led to Planck's law, the basics of quantum theory. Following the ideas of 
Planck on black body radiation, Einstein five years later (1905), quantized the electromag- 
netic field sj. Quantum thermodynamics is devoted to unraveling the intimate connection 
between the laws of thermodynamics and their quantum origin jjj J, [(3 
;he present study is aimed toward the quantum study of the third 



211 ] . In this tradition 



n 



23| . in particular quantifying the unattainability principle 24]: What is the scaling 



aw of thermodynamics 



of the cooling power J c of a refrigerator when the cold bath temperature approaches the 
absolute zero J c oc T a when T c — > 0. 



II. THE QUANTUM TRICKLE 



The minimum requirement for a quantum thermodynamical device is a system connected 
simultaneously to three reservoirs 25]. These baths are termed hot, cold and work reservoir 
as described in Fig. [1] A quantum description requires a representation of the dynamics 
working medium and the three heat reservoirs. A reduced description is employed in which 
the dynamics of the wor king medium is described by the Heisenberg equation for the operator 



O for open systems [26|, |27J : 



dt 



i ~ ~ 3d - 

-[H s , O] + — + C h (Q) + C c (0) + C W (Q) , 



(1) 



where H s is the system Hamiltonian and C g are the dissipative completely positive super- 
operators for each bath (g = h,c,w). A minimal Hamiltonian describing the essence of the 
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FIG. 1: The quantum trickle: A quantum heat pump designated by the Hamiltonian H s coupled 
to a work reservoir with temperature T w , a hot reservoir with temperature and a cold reservoir 
with temperature T c . The heat and work currents are indicated. In steady state Jh + J c + V = 0. 



quantum refrigerator is composed of three interacting oscillators: 
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frjjJint ( a^bc + ab^c^ 



(2) 



Hint represents an annihilation of excitations on the work and cold bath simultaneous with 
creating an excitation in the hot bath. In an open quantum system the superoperators L g 
represent a thermodynamic isothermal partition allowing heat flow from the bath to the 
system. Such a partition is equivalent to the weak coupling limit between the system and 
bath 



10|. The superoperators C g are derived from the Hamiltonian: 



H = H s + H h + H c + H w + H sh + H sc + H s 



(3) 



where H 9 are bath Hamiltonians and H S5 represent system bath coupling. Each of the 
oscillators is linearly coupled to a heat reservoir for example for the hot bath: H s /i = 
X s h(3iA. h + a^A/j) . Each reservoir individually should equilibrate the working medium to 
thermal equilibrium with the reservoir temperature. In general, the derivation of a thermo- 
dynamically consistent master equation is technically very difficult [28]. Typical problems 
are approximations that violate the laws of thermodynamics. We therefore require that 
the master equations fulfill the thermodynamical laws. Under steady state conditions of 
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operation they become: 

Jh + Jc + V = 



(4) 

2k. _ £ _ v > o , 



T^h -^c T w 



S3, 



where J7fe = (£fc(H)}. The first equality represents conservation of energy (first law) 
and the second inequality represents positive entropy production in the universe S M > 
(second law). For refrigeration T w > > T c . From the second law the scaling exponent 

a > 1 |n|. 

III. NOISE DRIVEN REFRIGERATOR 

Gaussian noise driven refrigerator. In the absorption refrigerator the noise source 
replaces the work bath and its contact hjj w c*c leading to: 

Hw = /WjVb + ab 1 ) =/(t)X, (5) 

where /(t) is the noise field. X = (a^b+ab^) is the generator of a swap operation between the 
two oscillators and is part of a set of SU (2) operators , Y = i(a^b — ab ), Z = ^a^a — b^b^ 
and the Casimir N = ^a^a + b^bj . 

We first study a Gaussian source of white noise characterized by zero mean (f(t)) = 
and delta time correlation (f(t)f(f)) = 2i]5(t — t'). The Heisenberg equation for a time 
independent operator O reduced to: 

jO = i{fl s ,d]+£ n (6) + £ h (d)+£ c (d) , (6) 
where H s = huh^R + hw c \)h. The noise dissipator for Gaussian noise is C n (0) = 



-t?[X,[X,6]] |29|. 



The next step is to derive the quantum Master equation of each reservoir. We assume 
that the reservoirs are uncorrelated and also uncorrelated with the driving noise. These 
conditions simplify the derivation of £h which become the standard energy relaxation terms 
driving oscillator u^a^a to thermal equilibrium with temperature and £ c drives oscillator 
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hubbub to equilibrium T c [27|. 



C h (6) = r h (iV A + l)(a+6a-i{ata,6} 
+ r h JV ft (a6a f - §{^,6} 

£ c (6) = r c (iv c + i) (^ob-ljb^o} 
+ r c iv c (botf - \ {bb\ 6} 



(7) 



In the absence of the stochastic driving field these equations drive oscillator a and b sepa- 
rately to thermal equilibrium provided that N h = (exp(|^) — and N c = (exp(|^r) — l) -1 . 
The kinetic coefficients I\/ c are determined from the baths density function [101 ] . 

The equations of motion are closed to the SU(2) set of operators. To derive the cooling 
current J c = (C c (hu c b^b)), we solve for stationary solutions of N and Z, obtaining: 

Cooling occurs for A" c > A^ =>■ ^ > jr. The coefficient of performance (COP) for the 
absorption chiller is defined by the relation COP = with the help of Eq. flH]) we obtain 
the Otto cycle COP |3o|: 



COP = < _?L_ . (9) 

A different viewpoint starts from the high temperature limit of the work bath T w based on 
the weak coupling limit in Eq. (0), ([2]), then: 



C w (6) = T W (N W + 1) (Wb6b f a- ij^abb^o} 
+ F W N W (Woa+b - \ jaa^b, 6 



(10) 



where N w = (exp(-^) — l) -1 . At finite temperature C w {0) does not lead to a close set of 
equations. But in the limit of T w — > oo it becomes equivalent to the Gaussian noise generator: 
C w (d) = —rj/2 ^[X, [X, 6]] + [Y, [Y, 6]] J, where rj = T W N W . This noise generator leads 
to the same current J c and COP as Eq. (jSJ) and (jUj). We conclude that Gaussian noise 
represents the singular bath limit equivalent to T w — > oo. As a result the entropy generated 
by the noise is zero. 

The solutions are consistent with the first and second laws of thermodynamics. The COP 
is restricted by the Carnot COP. For low temperatures the optimal cooling current can be 
approximated by J c ~ u c T c N c . Coupling to a thermal bosonic field such as electromagnetic 



or acoustic phonons field implies T c oc where d is the heat bath dimension. Optimizing 
the cooling current with respect to u c one obtains that the exponent a quantifying the third 
law J c oc T c a is given by a — d + 1 . 

Poisson noise driven refrigerator. Poisson white noise can be referred as a sequence 
of independent random pulses with exponential inter-arrival times. These impulses drive the 



coupling between the oscillators in contact with the hot and cold bath leading to 28|, l3l| : 

™ = (i/h)[H,6]-(i/h)\(t)[%6] 



where H is the total Hamiltonian including the baths. A is the rate of events and £ is the 
impulse strength averaged over a distribution P(£). Using the Hadamard lemma and the 
fact that the operators form a closed SU(2) algebra, we can separate the noise contribution 
to its unitary and dissipation parts, leading to the master equation, 

™ = (i/h)[H, 6] + (i/h)[n', 6] + c n (6) . (12) 

The unitary part is generated with the addition of the Hamiltonian H = heK. with the 
interaction 

e=~ [ dtP{i){2£/h-sin{2t/h)) . 



2 . 

This term can cause a direct heat leak from the hot to cold bath. The noise generator C n {p), 
can be reduced to the form C n (0) = — r)[K, [X, O]] , with a modified noise parameter: 

r] = j(l-J d£P(Z)cos(2Z/h) 

The Poisson noise generates an effective Hamiltonian which is composed of H and H , 
modifying the energy levels of the working medium. This new Hamiltonian structure has to 
be incorporated in the derivation of the master equation otherwise the second law will be 
violated. The first step is to rewrite the system Hamiltonian in its dressed form. A new set 
of bosonic operators is defined 

Ai = acos(#) + bsin(0) 

(13) 

A 2 = bcos(6') — asin(6 l ) , 
The dressed Hamiltonian is given by: 

h s = nn+klM + m]_a$a 2 , (14) 



where Q± = h — - ± \ j{ h — -) 2 + e 2 and cos 2 (6 ) ) = — — — Eq. (fT4]) impose the 



2 v 2 ' " y ' Q + -Q_ 

restriction, Q± > which can be translated to UhUJ c > e 2 . The master equation in the 

Heisenberg representation becomes: 

d6 



where 



(jf (i/n)[H s ,0]+A(0) + £ c (0) + £ n (0) , (15) 



c h (6) = 7iC 2 (Ai6a{ - HAiAj, 6}) 



+7 2 l c 2 (A t 1 6A 1 -±{A t 1 A 1 ,6}) 
+7 3 ft s 2 (A 2 6A$-±{A 2 A 2 ,6}) 
+7 4 h s 2 (A 2 6A 2 -±{A 2 A 2 ,6}) 



2 

£ c (6)= 7i c s 2 (A 1 6A t 1 -i{A 1 A t 1 ,6}) 



(16) 



+7 2 c s 2 (A t 1 6A 1 -i{A t 1 A l! 6}) 

+7 3 c c 2 (A 2 6A 2 -I{A 2 A 2 ,6}) 

+7^c 2 (A 2 6A 2 -±{a 2 a 2 ,6}) 

where s = sin(#) and c = 003(6*). And the noise generator: 



£ n (6) = -7 ? [W,[W,6]] , (17) 

where W = sin(20)Z + cos(20)X and a new set of operators which form an SU(2) algebra 
is defined: X = (AIA2 + A^Ai) , Y = i(A{A 2 - A 2 Ai) and Z = (A^Ai - A 2 A 2 ). The 
total number of excitations is accounted for by the operator N = (AjAi + A 2 A 2 ). 
The generalized heat transport coefficients become = 7 2 — 7i and C- = it ~ I3 f° r 



k = h,c. Applying the Kubo relation 32j, [33|: 7* = e m +^^k anc j _ e nn-p k ^k^ leads to 
the detailed balance relation: 



ji - 1 = Arfc 



In general £± is temperature independent and can be calculated specifically for different 
choices of spectral density of the baths. For electromagnetic or acoustic phonon field £± oc 
VL±. The heat currents Jh, J c and J n are calculated by solving the equation of motion for 
the operators at steady state and at the regime of low temperature, where cos 2 (8) ~ 1 and 
sin 2 (6) w 0. 



f = -§(ci + on - |(ci - e)z + (cw + e^f ) 
§ = -Kci + oz - jcci - on + (c^{ - e^) - 4^ 



(18) 
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Once the set of linear equations is solved the exact expression for the heat currents is 
extracted, J h = (c h (H s )\ J c = ^£ C (H S )^) and J n = //^(H^. For simplicity, the 
distribution of impulses in Eq. ( 1TT|) . is chosen as P(£) = S(£ — £o)- Then the effective noise 
parameter becomes: 



77= ~(l-cos(2£oA)) • (19) 



The energy shift is controlled by: 



e=~(2Z /h-sin(2Z /ti)) . (20) 
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FIG. 2: Entropy production = —Jk/Tk as a function of impulse £0 f° r the cold E c hot 
and the total entropy production S u = £/, + S c . T c = 10~ 3 , = 2., w c = T c , ujh = 10. A = uj c 
(X = Wc/10 (h = k = 1). 

Figure [2] shows a periodic structure of the heat current J7 C and the entropy production 
S c = —J c /T c with the impulse £o- The second law of thermodynamics is obtained by 
the balance of the large entropy generation on the hot bath compensating for the negative 
entropy generation of cooling the cold bath. The COP for the Poisson driven refrigerator 
is restricted by the Otto and Carnot COP: 

COP = 7^ L 7T<^—<^F- (21) 
The heat current J c is given by: 

N c _-N% 

Jc ~ (2r/)- 1 + (C)- 1 + (C)- 1 ' (22) 
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The scaling of the optimal cooling rate is now accounted for. The heat flow is maximized 
with respect to the impulse £o by maximizing rj Eq. f[T9"j) . which occurs for £ — n \ > 
(n = 1,2..). On the other hand the energy shift e 2 Eq. (120]) should to be minimized. The 
optimum is obtained when £o = f ■ The cooling power of the Poisson noise case Eq. (122]) is 
similar to the Gaussian one Eq. (jHJ). In the Poisson case also the noise driving parameter rj 
is restricted by u c . This is because e is restricted by > and therefore A is restricted to 
scale with u c . In total when T c — > 0, J c oc T c d+1 . 

The optimal scaling relation J c oc T c a of the autonomous absorption refrigerators should 
be compared to the scaling of the discrete four stroke Otto refrigerators 34j. In the driven 
discrete case the scaling depends on the external control scheduling function on the expan- 
sion stroke. For a scheduling function determined by a constant frictionless nonadiabatic 
parameter the optimal cooling rate scaled with a = 2. Faster frictionless scheduling pro- 
cedures were found based on a bang-bang type optimal control solutions. These solutions 
led to a scaling of a = 3/2 when positive frequencies were employed and J c oc — T c /logT c 



35 



The drawback of the externally 



when negative imaginary frequencies were allowed 
driven refrigerators is that their analysis is complex. The optimal scaling assumes that the 
heat conductivity V 3> u c , and that noise in the controls does not influence the scaling. For 
this reason an analysis based on the autonomous refrigerators is superior. 
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